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A black hole mass threshold from non-singular quantum gravitational collapse 
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Quantum gravity is expected to remove the classical singularity that arises as the end-state of 
gravitational collapse. To investigate this, we work with a toy model of a collapsing homogeneous 
scalar field. We show that non-perturbative semi-classical effects of Loop Quantum Gravity cause a 
bounce and remove the black hole singularity. Furthermore, we find a critical threshold scale, below 
which no horizon forms - quantum gravity may exclude very small astrophysical black holes. 
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Singularity formation during gravitational collapse sig¬ 
nals the breakdown of classical general relativity. In a 
more complete theory of quantum gravity the singular¬ 
ity should be removed. However, a satisfactory quantum 
gravity theory has yet to be developed. In addition, the 
dynamics of general collapse is very complicated. Thus 
we can only expect to make partial progress in tackling 
the problem, using a candidate for quantum gravity and 
a collapse model that is simple enough to be tractable. 

A non-perturbative approach to quantizing gravity is 
loop quantum gravity or quantum geometry lj, which 
gives rise to a discrete spatial structure [5] and whose 
successes include prediction of black hole entropy Jj. Ap¬ 
plied to the early universe, loop quantum effects can re¬ 
move the big bang singularity {4j. A natural question 
is: do these effects also remove the black hole singular¬ 
ity as the end-state of collapse? Techniques to handle 
inhomogeneous systems are under development and give 
promising indications @|, but they do not easily reveal 
the physical picture. We thus consider a simple toy model 
of a collapsing homogeneous scalar field. Classically, this 
model always produces a black hole, but we show that 
loop quantum effects change this situation dramatically. 

Since we do not yet know semi-classical non- 
perturbative effects in inhomogeneous cases, we are un¬ 
able to perform our analysis in the general case. However, 
when we split the system into a homogeneous star inte¬ 
rior and an inhomogeneous outside region, known quan¬ 
tum effects in the interior can be carried into the ex¬ 
terior indirectly through matching conditions. The col¬ 
lapsing homogeneous scalar field cannot be matched to 
a Schwarzschild exterior because the pressure does not 
vanish at the boundary. But in any case, we expect that 
quantum effects will include small non-stationary correc¬ 
tions and thus use a non-stationary spherically symmetric 
exterior. The generalized Vaidya metric provides a rea¬ 
sonable starting point. It is sufficiently general to allow 
for a broad range of behavior, including radiative effects. 

Our analysis is based on effective equations for the in¬ 
terior which have been established in the cosmological 


setting. Fundamentally, the evolution is described by a 
wave function subject to a difference equation, and effec¬ 
tive equations describe the motion of semiclassical wave 
packets Q. As long as one stays in semiclassical regimes, 
which can e.g. be checked using the size of curvature, one 
gets reliable expectations for the quantum situation. 

We first review the classical collapse and the inevitabil¬ 
ity of a black hole singularity covered by a horizon, for 
any initial mass. The isotropic interior metric is @] 

ds 2 = —dt 2 + a(t) 2 (l + r 2 /4) -2 [dr 2 + r 2 dfl 2 ] , (1) 

and the massless scalar field <fi(t) has pressure and energy 
density p = p = ^<j) 2 . The Friedmann equation is 

a 2 /a 2 = 47t^ 2 c/> 2 /3 — 1/a 2 . (2) 

The Klein-Gordon equation, acf> + 3d <f> = 0, has solution 
4 > = L/a 3 , (3) 


where L is a length scale associated with the maximal 
size of the collapse region, since 0 implies 

a < a m = (47t/3) 1/4 \JTpL. (4) 


At the singularity a —> 0, we have (f>,p —> oo. The solu¬ 
tion of the Friedmann equation is 


t to — a m 
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where ao(< a m ) gives the initial size of the collapse region 
at time to- The singularity a = 0 is covered by a horizon 
(see below), and reached in finite proper time for any ag: 
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where F is an elliptic integral of the first kind. 
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We now add non-perturbative modifications to the dy¬ 
namics, motivated by loop quantum gravity The 
quantization introduces a fundamental length scale 

4 = 0.28 VH, (7) 


where j(> 1) is a half-integer that is freely specifiable. 
For a < 4, the dynamics is increasingly different from 
general relativity. For a < £ p , the continuum approxima¬ 
tion to the spacetime geometry begins to break down, 
and the fully quantum gravity regime is reached. In 
the intermediate regime £ p < a < 4, loop quantum ef¬ 
fects may be treated semi-classically, i.e., the spacetime 
metric behaves classically, while the dynamics acquires 
non-perturbative modifications to general relativity [f|. 
The non-perturbative semi-classical regime exists pro¬ 
vided 4 £ p , i.e., for j 1. 

The key feature of the loop quantization scheme is 
the prediction that the geometrical density, 1/a 3 , does 
not diverge as a —> 0, but remains finite. The expec¬ 
tation values of the density operator are approximated 
by dj(a) = D(a) a -3 , where the loop quantum correction 
factor is i,9] 


D(a) = (8/77) 6 q 3,2 [l\{q + 1) 11/4 - I q- 1| 11/4 
- 11 q[(q + 1) 7/4 - Sgn (q - 1 )|q - 1| 7 / 4 ] (8) 


with q = a 2 /£ 2 . In the classical limit we recover the ex¬ 
pected behavior of the density, while the quantum regime 
shows a radical departure from classical behavior: 

a > f» : D « 1 , a<4 :f~(12/7) 6 (a/4) 15 . (9) 


Then dj remains finite as a —> 0, unlike in conventional 
quantum cosmology, thus evading the problem of the big- 
bang singularity in a closed model jlOj . Intuitively, one 
can think of the modified behavior as meaning that grav¬ 
ity, which is classically always attractive, becomes re¬ 
pulsive at small scales when quantized. This effect can 
produce a bounce where classically there would be a sin¬ 
gularity, and can also provide a new mechanism for high- 
energy inflationary acceleration EH In the semi-classical 
regime (where the spectrum can be treated as continu¬ 
ous), dj has a smooth transition from classical to quan¬ 
tum behavior, varying from a -3 to a 12 . We emphasize 
that this is but one possibility for a bounce which we use 
for concreteness, while bounces in general appear more 
generically in loop cosmology E3 

In loop cosmology the Hamiltonian of a scalar field in 
a closed universe is 


n = a 3 v(t) + dj pH 2, Pj, = d-\ 


( 10 ) 


where Pj, is the momentum canonically conj ugat e to <j>. 
This leads to a modified Friedmann equation 111 1131. 
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FIG. 1: The scale factor a(t) of the collapsing interior, 
for classical (dashed) and semi-classical quantum dynamics 
(solid). 


and a modified Klein-Gordon equation El 

c/+3 a _1 a (1 —a) 4>+DV{(j)) = 0 , a = aD/(3dD). (12) 

For a <C 4, we have a —♦ 5, whereas classically D = 1 
and hence a = 0. Thus in the semi-classical regime, 
0 < a < 5. 

For a massless scalar field, V = 0, the solution of 
Eq. <11211 . generalizing Eq. <0, is 

</> = Ldj(a) , (13) 

so that Pj, = L = const. Then the Friedmann equation 
becomes 

d 2 + 1 = D(a)(a m /a) 4 . (14) 

The energy density and pressure are modified as p = 
<j) 2 /2D , p = <j) 2 (1 — a) /2 D , so that 

w=p/p=l — a. (15) 

(The modified p and p satisfy the usual conservation 
equation if (f> satisfies the modified Klein-Gordon equa¬ 
tion.) Since a varies from 0 to 5 as a decreases, the 
c[> term in Eq. o. which classically behaves as anti- 
frictional during collapse, starts to behave as frictional 
when a > 1. Thus, contrary to classical behavior, where 
(f> increases as a decreases, in the semi-classical regime 
the scalar field starts slowing down with collapse. In 
fact at a = 2 the magnitude of the frictional term be¬ 
comes exactly equal to the classical anti-frictional term. 
Thereafter at smaller values of the scale factor the term 
becomes increasingly frictional and the collapse further 
slows down, and may turn around. 

The point where a = 2 is also the point beyond 
which the null energy condition is violated: w < — 1, 
by Eq. JG3- Violations of the null energy condition by 
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quantum gravity effects are to be expected, and in loop 
quantum gravity this occurs for a > 2, when the scalar 
field effectively behaves as a “phantom” field. 

In order to see qualitatively how the non-perturbative 
frictional quantum effects remove the classical singular¬ 
ity, we assume that, over a small interval of scale factor, 
we can take a ~ constant, so that D « ID* (a/I*) 3a , 
where ID* is a dimensionless constant. By Eq. 

0 « LID*C 3a « 3(a_1) , (16) 

which shows how the kinetic energy decreases with de¬ 
creasing a when a > 1, contrary to the classical case. 
The modified Friedmann equation m gives 

d 2 « «C 3q ^*) a 3 “ -4 - 1 • (17) 

In general relativity, where a = 0 and ID* = 1, this shows 
that for a < a m , there is no turning point in a, i.e., a ^ 0. 
With loop quantum effects, for a > |, the equation 

a(t c ) = 0 has a solution, a c ~ (t 3a / 1 ^ 3a 4 ' ) -C 

a m . Thus the collapse leads to a bounce and singular¬ 
ity avoidance. The numerical integration of the mod¬ 
ified Friedmann and Klein-Gordon equations confirms 
the qualitative analysis, and the results are illustrated 
in Fig. 0 As is clear from the figure, the classical curve 
(dashed line) hits the singularity in finite time, whereas 
the quantum-corrected curve bounces and avoids the sin¬ 
gularity. The key question is whether a horizon forms in 
the quantum-corrected collapse. 

The formation or avoidance of the singularity a = 0 
is independent of the matching to the exterior. But 
in order to understand horizon formation in the semi- 
classical quantum case, we need to impose the match¬ 
ing conditions. Since the pressure is nonzero at the 
boundary, given in comoving coordinates by r = R = 
constant, the interior cannot be directly matched to a 
static Schwarzschild exterior. However we can match to 
an intermediate non-stationary region - for example, a 
generalized Vaidya region 0, 

ds 2 = — [1 — 2 M(v, x)/x] dv 2 + 2 dvd% + x 2 dfl 2 . (18) 

The usual Vaidya mass M /££ is generalized so that 
dM/dx may be nonzero. The total mass measured by 
an asymptotic observer is to = tom + Wf where tom 
is the total mass in the generalized Vaidya region, and 
= f pdV the interior mass. By Eqs. 0, 0 and 
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(19) 


Since we do not specify the matter content in the exterior, 
and do not know the modified field equations there, we 
cannot determine M(v , x) and thus tom- However, as we 
discuss below, we can still draw qualitative conclusions 
about the behavior of horizons close to the matter shells. 



FIG. 2: The speed of collapse, |a|, against the scale fac¬ 
tor a, for the evolution shown in Fig. 0 up to the bounce. 
The dashed curve is for classical dynamics and semi-classical 
quantum dynamics gives the solid curve. The horizontal dot¬ 
ted lines correspond to different values of R in Eq. 1241 : for 
the upper line there is no horizon in the quantum-corrected 
case, the middle line corresponds to the threshold for a hori¬ 
zon, the lower line to the case of an inner and outer horizon. 


Matching the 1st and 2nd fundamental forms we obtain 

X (v) = Ra(t)/(1 + H 2 /4), (20) 

dv/dt = (1 + H 2 /4)/(l — H 2 /4 — Ra) , (21) 

2 M = aR 3 {d 2 + l)/(l + H 2 /4) 3 , (22) 

-M, v = x, w + (1-2M/x-3x,„)(M/x-M )X )(23) 

The exterior region can contain trapped surfaces when 
the condition 2 M(v,x) = X is satisfied. Evaluating this 
at the matching surface, using Eqs. m and (E3, gives 

|a| = H -1 (l — H 2 /4). (24) 

When this value is reached, a dynamical horizon 0 in¬ 
tersects the matching surface. This always occurs clas¬ 
sically since during the collapse |a| varies from zero to 
infinity. With the modified dynamics, however, |d| is 
bounded throughout the evolution, so that it depends on 
initial values whether or not a horizon forms (Fig. 0. 
Moreover, after the bounce, a grows again, so that the 
condition can be satisfied a second time. This results in 
a picture where the bounce, replacing the classical sin¬ 
gularity, may be shrouded by an evaporating dynamical 
horizon outside, as shown in Fig. 0 There will be a sec¬ 
ond point where the horizon condition is satisfied since 
|a| decreases between the peak of dj(a) and the bounce. 

When it intersects the matching surface, the horizon is 
always null, as follows from Eq. (E3>. Its later behavior 
depends on the details of the outer region, which can not 
be determined here. Nevertheless, one can expect that 
both horizons will become timelike and evaporate. Hori¬ 
zon evaporation in this model does not only come from 
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FIG. 3: Eddington-Finkelstein diagram of the collapse, with 
boundary E. Dotted lines show constant v (outside) and con¬ 
stant t (inside). Quantum modifications imply a bounce of 
the collapsing field, which for large enough mass is covered 
by an inner and outer evaporating horizon (dashed). A sin¬ 
gle matching suffices only until the inner horizon disappears. 
The dot-dash curves correspond to the subsequent evolution 
which is not determined in our model. 


Hawking radiation, which may be included effectively in 
the outside matter content, but also from violations of 
energy conditions around the bounce, which may lead to 
effective outgoing negative energy. 

The model is not able to specify the future of the sys¬ 
tem after it re-emerges out of the horizon. Equation m 
shows that dv/dt diverges if and only if a > 0 and the 
matching surface becomes trapped. Thus, we can de¬ 
scribe the collapse with a single matching until a horizon 
disappears, at which point the interior t ceases to be 
a good coordinate. One has to continue with a second 
matched region to analyze the future of the system, but 
this is beyond the scope of our model. 

The qualitative picture that emerges from our toy 
model is thus the following: 

• We do obtain black holes, i.e., dynamical horizons, for 
large masses, but they contain a bounce of the infalling 
matter rather than a singularity. For large mass, vio¬ 
lations of energy conditions are initially small and the 
evaporation takes a long time, so that there are only small 
deviations from classical results. 

• For small enough mass however, black holes do not 
form; horizons do not develop during collapse and the 
bounce is uncovered. The critical threshold scale for 
horizon formation is given by the turning point in the 
|d| curve. By Eqs. m and GJ, the critical scale is 

(25) 


The corresponding threshold mass is m cr it = Tom + 
m<t>(a C iit), but we are unable to compute this mass be¬ 
cause the exterior dynamics remains undetermined. 

Our estimates may be strongly influenced by the sim¬ 
plifications, in particular a homogeneous interior, we are 
forced to impose on the problem. However, the qualita¬ 
tive features should be robust, and can provide guidance 
for further more general analysis. In particular, they 
mean that there could be lower bounds on the masses 
of black holes that form by gravitational collapse. This 
could rule out primordial black holes below the thresh¬ 
old mass, and thus modify estimates of Hawking radi¬ 
ation effects from very small black holes. More specu¬ 
lative is an extension to highly non-spherical situations 
such as particle collisions. If LQG effects can in fu¬ 
ture be shown to encode some of the non-perturbative 
aspects of string theory, then our results may have im¬ 
plications for the production of black holes in colliders, 
as predicted in brane-world gravity E3- The black hole 
horizon threshold would be a multiple not of £ p , but of 
the higher- dimensional Planck scale, which could be as 
low as O(TeV). This would mean that higher collision 
energies are needed for horizon formation, so that black 
hole production could be significantly reduced. 
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